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We measure directed flow (wi) for charged particles in Au+Au and Cu+Cu collisions at ^/sj^ ~ 
200 GeV and 62.4 GeV, as a function of pseudorapidity {rj), transverse momentum (pt) and collision 
centrality, based on data from the STAR experiment. We find that the directed flow depends on the 
incident energy but, contrary to all available model implementations, not on the size of the colliding 
system at a given centrality. We extend the validity of the limiting fragmentation concept to vi 



in different collision systems, and investigate possible explanations for the observed sign change in 
vi{pt). 

PACS numbers: 25.75.Ld 



The heavy ion program at the Rclativistic Heavy Ion 
Collider (RHIC) seeks to understand the nature and dy- 
namics of strongly-interacting matter under extreme con- 
ditions. It is widely expected that in collisions at RHIC, a 
new partonic phase of matter is created. sQGP, strongly 
interacting Quark Gluon Plasma P|. In particular, its 
bulk nature is revealed in strong elliptic flow, which in 
central collisions approaches the predictions of ideal hy- 
drodynamics, assuming system thermalization on an ex- 
tremely short timescale (~ 0.5 fm/c) [31 • However, the 
mechanism behind such rapid thermalization remains far 
from clear and is under active theoretical study [3|, Ij, l5| . 
This may be related to another novel phenomenon that 
could be relevant at RHIC — saturation of the gluon 
distribution — which characterizes the nuclear parton 
distribution prior to collision [61] . Various theoretical ap- 
proaches to connect collision geometry, saturated gluon 
distributions, and the onset of bulk collective behavior 
are being explored [3]; more experimental input would 
guide these efforts. 

Directed flow refers to collective sidewards deflection 
of particles and is characterized by a first-order har- 
monic (vi) of the Fourier expansion of particle's az- 
imuthal distribution w.r.t. the reaction plane [7|. At 
large r] (in the fragmentation region) the directed flow is 
believed to be generated during the nuclear passage time 
(2R/j ~ 0.1 fm/c) lalsl. It therefore probes the onset of 
bulk collective dynamics during thermalization, provid- 
ing valuable experimental guidance to models of the pre- 
equilibrium stage. In this Letter, we present multiple- 
differential measurements of vi for Au-I-Au and Cu-I-Cu 
collisions at ^/sNN =200 and 62.4 GeV as a function of 
?7, pt, and collision centrality. Here, we report an in- 
triguing new universal scaling of the phenomenon with 
collision centrality. Existing implementations of Boltz- 
mann/cascade and hydrodynamic models are unable to 
explain the measured trends. 

At RHIC energies, it is a challenge to measure vi accu- 
rately due to the relatively small signal and a potentially 
large systematic error arising from non-flow (azimuthal 
correlations not related to the reaction plane orienta- 
tion). In this work, the reaction plane was determined 
from the sideward deflection of spectator neutrons @, [lfl| 
measured in the Shower Maximum Detectors (SMD) 
of the Zero Degree Calorimeters (ZDC) [ul [12]. The 
vi based on this quantity, denoted i;i{ZDC-SMD} [ll| . 
should have minimal contribution from non-flow effects 
due to the large rj gap between the spectator neutrons 
used to establish the reaction plane and the 77 region 
where the measurements were performed. 

Charged particle tracks were reconstructed in STAR's 
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FIG. 1: (color online) Charged particle vi{'i]) for three central- 
ities in Au-f An collisions at 200 GeV. The arrows indicate the 
algebraic sign of vi for spectator neutrons, and their positions 
on the 77 axis correspond to beam rapidity. The inset shows 
the mid-77 region in more detail. The error bars are statisti- 
cal, and the shaded bands show systematic errors. PHOBOS 
results 1181 are also shown for mid-central collisions. 
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with pscudora- 
< 4.0, respec- 



main TPC [IJ and forward TPCs 
pidity coverage \t]\ < 1.3 and 2.5 < ]77] 
tively. The centrality definition (in which zero repre- 
sents the most central collisions) and track quality cuts 
are the same as in Ref. 15|. This study is based on 



Au+Au samples of eight million events at 200 GeV, five 
million at 62.4 GeV, and Cu-|-Cu samples of twelve mil- 
lion events at 200 GeV, and eight million at 62.4 GeV. 
All were obtained with a minimum-bias trigger. System- 
atic uncertainties on vi measurements are estimated to 
be within 10% for the 77 range studied. This limit is based 
on comparisons of wi{ZDC-SMD} and independent anal- 
ysis methods [ll|, ll5| , and we also make use of forward- 
backward symmetry to constrain estimates of systematic 
errors. Non-flow is not the dominant source of system- 
atic uncertaint y. M ore details about these errors can be 
found in Refs.|ll|,[li, 

The resolution [7|] of the first-order event plane recon- 
structed using the ZDC-SMDs is a crucial quantity for 
this analysis. The magnitude of the event plane resolu- 
tion, defined as {cosI^Pep — ^rp)) [2|, increases with the 
spectator vi and the number of neutrons per event de- 
tected by the ZDC-SMDs. The ZDC size is optimized for 
200 GeV, and its acceptance for spectator neutrons de- 
creases at lower energies due to spectator neutrons being 
emitted within a cone whose apex angle increases with 
the inverse of the beam momentum. For the 30 — 60% 
most central collisions, resolutions for 200 GeV Au+Au 
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FIG. 2: (color online) Charged particle vi{pt) in 200 GeV 
Au+Au for three centralities. The dashed and dotted curves 
are hydrodynamic calculations for the labeled rapidities at 
impact parameter 6.8 fm (15 — 25% most central collisions). 
See the text for an explanation of the solid curve. The plotted 
error bars are statistical, and systematic errors (see Fig.[T)) are 
within 10%. 



and Cu+Cu, and for 62.4 GeV Au+Au and Cu+Cu are 
about 0.4, 0.15, 0.15 and 0.04, respectively (more details 
are provided in Table 1 of Ref. ^.) The 30 - 60% cen- 
trality interval is the only region where the ZDC-SMD 
event-plane resolution can be reliably determined for all 
four systems. 

The charged-particle vi{t]) is shown in Fig. [1] for 
Au+Au at ^/s^ = 200 GeV in three centrahties. The 
inset shows, on expanded scales, the mid- 77 region mea- 
sured by the main TPC, where vi is resolvable below 
the 0.1% level. Within the studied 77 range, the sign of 
charged particle f 1 is opposite to that of the spectators, 
and the vi magnitude increases from central to peripheral 
collisions. For 0-5% centrality, the slope dvx/drj changes 
sign above the middle of the FTPC pseudorapidity accep- 
tance, and our results agree with t he p attern reported by 
PHOBOS over a broader -q range [lEEll- 

The ratio {px) / (pt) is shown in Fig.[T]for the most cen- 
tral data (0 to 5%), in comparison to vi. Here, Px refers 
to the in-plane component of a track's transverse momen- 
tum, a quantity commonly used prior to the 1990s [l0|. 
As elaborated below, there is interest in the behavior of 
both vi and (px) when vi{pt) changes sign. 

To further examine vi, the 200 GeV Au+Au data are 
divided into bins of pt (Fig. U). The upper and lower 
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FIG. 3: (color online) Charged particle vi{ri) for mid-central 
(30 - 60%) Au+Au and Cu+Cu at 200 GeV and 62.4 GeV. 
The solid and dashed curves are odd-order polynomial fits, 
to guide the eye and demonstrate the forward-backward sym- 
metry of the data. The wider shaded bands are from AMPT 
for the same conditions as the data. For clarity, 200 (62.4) 
GeV calculations are shown only at negative (positive) 77. The 
plotted error bars are statistical, and systematic errors (see 
Figs. [1] and [5} are within 10%. 



panels show results from the main TPC and the FTPCs, 
respectively. In the main TPC, wi(pt) crosses zero at 
I < Pt < 2 GcV/c for central and mid-central collisions. 
A zero-crossing behavior in vi{pt) is necessarily exhibited 
by a hydrodynamic calculation in which (px), presumably 
imparted during the passing time of the initial-state nu- 
clei, has been neglected and set equal to zero [19| . Due to 
the poor momentum resolution of the FTPCs at higher 
Pt, wc cannot test the zero crossing at forward ry. It is 
noteworthy that the observed lj>x), presented in Fig. [U 
is far from negligible, which contradicts the assumptions 
used in the hydrodynamic calculations. 

The observed wi(pt) dependence can be explained by 
assuming that pions and baryons flow with opposite 
sign, coupled with the measured baryon enhancement at 
higher pt [20|. For example, taking linear functions [21j 
for pion and baryon vi{pt), we obtain a satisfactory de- 
scription of our data (see the solid curve in Fig. [2| with 
pion vi slopes, dvi/dpt = -0.18 ± 0.02, -0.34 ± 0.02 
and -0.52 ± 0.04, and baryon vi slopes 0.56 ± 0.12, 
0.86 + 0.10 and 1.02 + 0.12 for centralities 0-5%, 5-40% 
and 40 — 80%, respectively. Note that the opposite vi 
slope for pions and protons, with the magnitude of pro- 
ton slopes being larger, in this case is consistent with 
calculations [221 where the "wiggle" rapidity dependence 
of identified particles has been predicted to result from 
the interplay of stopping and radial flow. Currently, we 
arc unable to test the wiggle effect in vi{y) with identi- 
fied particles due to limited statistics and limited particle 
identification. 



To study the energy and system size dependence of vi , 
Fig. [3] shows Cu+Cu data compared to Au+Au in the 
centrahty range 30 - 60% for both 200 and 62.4 GeV. 
There is a clear trend for vi (rj) to decrease with increasing 
beam energy for both Au+Au and Cu+Cu. In the stud- 
ied pseudorapidity and centrahty range, vi{ri) is, within 
errors, independent of the system size at each beam en- 
ergy, despite the three-to-one mass ratio between gold 
and copper. This remarkable feature holds for almost all 
centrahty bins studied, as shown in Fig. [H and persists 
even near mid-ry (as shown in the upper panel), where 
elliptic flow {V2) of charged particles in Cu+Cu is con- 
siderably lower than in Au+Au [23[. Unlike W2/e, the 
ratio of the elliptic flow to the system initial eccentricity, 
which scales with the particle density in the transverse 



plane {l/S)dNch/dy 24| (also interpreted to be the mid- 
rapidity area density 25 1 or the system length [26|), vi{ri) 
at a given centrahty is found to be independent of the sys- 
tem size, and varies only with the incident energy. The 
different scalings for V2/e and vi might arise from the way 
in which they are developed: to produce V2, many mo- 
mentum exchanges among particles must occur (and the 
number of momentum exchanges is related to the partic- 
ipant density and the dimensions of the system), while to 
produce vi , an important feature of the collision process 
is that different rapidity losses need to occur (related to 
the incident energy) for particles at different distances 
from the center of the participant zone |22||. 

The hybrid transport model AMPT [27[ lies consis- 
tently below the measured data, as evident from Fig. [3l 
STAR'S prior vi study ^ in Au+Au at 62 GeV also 
showed this trend for AMPT and other transport mod- 
els. It is noteworthy that AMPT does not exhibit the 
observed pattern of system-size independence. UrQMD 
[28| (not shown here) is similar to AMPT in exhibiting a 
significant change in vi between Au+Au and Cu+Cu. 

Further scaling behavior is seen by transforming the 
data presented in Fig. [3] into the projectile frame (see 
Fig. [5]), where zero on the horizontal axis corresponds 
to the beam rapidity, yboam, for each of the collision 
energies. Within three units from j/bcam, most data 
points lie on a universal curve for vi versus ij — j/boam- 
This incident-energy scaling of directed flow has previ- 
ously been reported for Au+Au 



11, hi 



and it is now 
evident that the limiting fragmentation hypothesis [29| 
holds even for much lighter collision systems like Cu+Cu. 
AMPT adheres less closely to limiting fragmentation for 
Cu+Cu. Note that the quantity r/ — j/bcam introduces 
some uncertainty due to the use of rj instead of rapidity; 
the latter requires particle identification. The system- 
size independence at a given fractional cross section and 
longitudinal scaling of scaled multiplicity distributions, 
dNch/di]/{Npart/'2), have been previously reported by 
the PHOBOS Collaboration [s^. 

In summary, we have presented measurements of 
charged-particle directed flow as a function oi pt, 77 and 
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FIG. 4: (color online) Charged particle vi versus centrality, 
for Au+Au and Cu+Cu at 200 GeV and 62.4 GeV. The upper 
(lower) panels show results from the main TPC (FTPC). The 
plotted error bars are statistical, and systematic errors (see 
Figs. [1] and [5} are within 10%. 
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FIG. 5: (color online) Charged particle vi versus 77 — yheam, 
for 30 - 60% Au+Au and Cu+Cu at 200 and 62.4 GeV. The 
plotted error bars are statistical, and the shaded bars show 
systematic errors. 



centrality in Au+Au and Cu+Cu collisions at ^s„„ — 
200 and 62.4 GeV. The observed trend of decreasing vi 
with increasing beam energy agrees with models. The 
lack of system-size dependence in vi for Au+Au and 
Cu+Cu is quite remarkable and is a feature not observed 
or predicted by any existing model implementation. The 



presented 77-dependence of vi provides further support 
for limiting fragmentation scaling by extending its ap- 
plicability to Cu+Cu. The observed pt-dependence of 
directed flow motivates further theoretical investigations 
and experimental measurements with identified particles. 
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